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crack  between  two  materials  in  a  composite  is  oscillatory  when  the  two 
materials  are  mismatched  .  This  leads  to  the  physically  unacceptable 
interpenetration  of  the  crack  surfaces.  However  it  is  saved  from  a 
higher  order  stress  singularity  when  the  effects  of  heat  flow  is  included. 
On  the  other  hand,  when  the  two  materials  are  "not  mismatched",  the  inter¬ 
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Final  Report  on 
AFOSR  Contract  F49620 


Anisotropic  Elastic  Materials  and  Composites  With 
The  Effects  of  Temperature  and  Piezoelectricity 


The  objectives  of  this  project  were  to  extend  the  Stroh  formalism  for 
two-dimensional  deformations  of  anisotropic  elastic  bodies  to  include  heat 
flow  and  piezoelectricity,  and  to  apply  the  theory  to  study  the  effects  of 
temperature  and  piezoelectricity  in  composites.  The  basic  Stroh  formalism 
is  not  new.  The  foundation  of  the  formalism  was  laid  down  by  Stroh  in  his 
pioneering  papers  in  1958  and  1962.  Since  then  material  scientists,  physi¬ 
cists,  and  applied  mathematicians  have  recognized  its  superior  features  and 
have  further  developed  the  formalism.  For  unexplainable  reasons  the 
engineering  community  has  refused  to  employ  Stroh  formalism  and  have 
steadfastly  employed  the  traditional  Lekhnitskii  formalism.  Only  in  recent 
years  did  a  few  young  researchers  in  the  applied  mechanics  community 
begin  to  employ  the  Stroh  formalism.  The  Stroh  formalism  is  mathemati¬ 
cally  elegant  and  technically  powerful.  It  is  not  difficult  to  understand.  In 
fact  it  is  easier  to  learn  and  simpler  to  use  than  the  Lekhnitskii  formalism. 
In  view  of  this,  it  is  even  more  puzzling  why  the  engineering  community 
has  been  so  slow  in  catching  up  with  the  Stroh  formalism.  Researchers  in 
material  sciences,  physics,  and  applied  mathematics  hardly  use  the  Lekhnit¬ 
skii  formalism. 
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In  the  first  year  of  the  project  we  extended  the  Stroh  formalism  to 
include  heat  flow  in  general  anisotropic  elastic  materials  [1].  We  then 
studied  the  stress  singularities  at  an  interface  crack  due  to  heat  flow  in  a 
bimaterial  that  consists  of  two  dissimilar  general  anisotropic  materials.  If 
the  stress  singularities  at  the  crack  tip  are  written  in  the  form  of  rs,  there 
are  three  stress  singularities 

r  2,  r  2+<r,  and  r  2  'y  (1) 

when  the  heat  effects  are  ignored.  The  materials  in  the  composite  are 
mismatched  when  y  0  so  that  the  second  and  the  third  8  in  (1)  are 
complex.  With  the  inclusion  of  heat  flow,  it  is  shown  that  there  is  a 
possibility  of  a  fourth  stronger  singularity 

r  2 (In r) .  (2) 

This  is  a  stronger  singularity  than  the  ones  shown  in  (1),  indicating  that 
heat  flow  increases  the  potential  for  initiation  of  fracture  in  composites. 
The  conditions  for  the  existence  of  the  stronger  singularity  (2)  are 
presented  explicitly. 

The  conditions  for  the  existence  of  the  stronger  singularity  derived  in 
[1]  are  applicable  to  general  anisotropic  materials.  It  is  difficult  to 
interpret  the  conditions  physically.  In  order  to  have  a  quantitative  result, 
monoclinic  materials  with  the  symmetry  planes  at  x3=0  are  considered  [2]. 
All  quantities  involved  in  the  solutions  and  in  the  conditions  for  the 
existence  of  the  stronger  singularity  now  have  an  explicit  real  form 
expression.  An  unexpected  result  is  that  the  stronger  singularity  does  not 
exist  if  the  two  materials  in  the  bimaterial  are  mismatched.  The  is  a  rather 
puzzling  result.  One  would  have  expected  that  the  reverse  is  true.  It  is 
known  that  a  mismatched  bimaterial  induces  an  oscillatory  displacement 
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near  the  crack  tip,  resulting  in  the  physically  unacceptable  interpenetration 
of  the  two  materials  at  the  interface  crack.  It  is,  however,  saved  from  a 
stronger  stress  singularity  when  the  effects  of  heat  flow  is  included.  On 
the  other  hand,  while  a  non-mismatched  bimaterial  is  immune  to  the 
interpenetration  of  the  two  materials  at  the  interface  crack,  it  may  admit  a 
stronger  stress  singularity  when  the  heat  flow  in  considered. 

For  isotropic  bimaterials  the  two  materials  are  not  mismatched  if 

1  -2v  1  -  2  v' 


(3) 


where  n  is  the  shear  modulus  and  v  is  the  Poisson’s  ratio.  The  prime  refers 
to  the  second  material.  The  condition  for  existence  of  the  stronger 
singularity  is 


a(l+v)  t«(l+v') 

- * - ; —  (4) 

K  K 

where  a  and  k  are  the  thermal  expansion  coefficient  and  the  heat  conduc¬ 
tivity,  respectively.  Thus  when  (3)  and  (4)  are  satisfied,  there  exists  a 
stronger  stress  singularity  of  the  form  (2)  in  an  isotropic  bimaterial.  If  the 
material  is  incompressible,  i.e.,  if  v=l/2,  (3)  is  automatically  satisfied  and 
(4)  leads  to 


cck'^cc'k.  (5) 

In  the  second  and  third  year  we  looked  at  the  extension  of  the  Stroh 
formalism  to  piezoelectric  materials  and  its  applications.  The  sextic  formal¬ 
ism  of  Stroh  becomes  an  octet  formalism.  Most  of  the  identities  for 
anisotropic  materials  can  be  extended  to  piezoelectric  materials  except  that 
certain  matrices  are  no  longer  positive  definite.  However,  they  can  be 
shown  to  be  nonsingular,  a  very  important  and  crucial  property  in  the 
analysis  of  piezoelectric  materials  [3]. 
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One  of  the  basic  problems  is  that  of  the  Green's  function  for  the  infinite 
space  subject  to  a  prescribed  singularity  at  the  origin  of  the  coordinate 
system.  The  usefulness  of  Green's  function  in  solving  practical  problems 
with  an  arbitrarily  prescribed  boundary  conditions  on  an  irregular  shape  of 
boundary  is  well  recognized.  We  presented  in  [3]  the  Green's  function  in 
the  infinite  space  of  homogeneous  piezoelectric  material  subjected  to  a  line 
force,  line  charge,  and  a  line  dislocation.  To  make  it  more  general,  we 
consider  the  infinite  space  to  consist  of  an  arbitrary  number  of  wedges  of 
different  wedge  angles  and  different  piezoelectric  materials.  The  Green’s 
function  due  to  a  line  force,  a  line  charge,  and  a  line  dislocation  at  the 
center  of  this  composite  space  consisting  of  an  arbitrary  number  of  wedges, 
has  a  surprisingly  simple  expression.  Moreover,  the  solution  can  be 
expressed  in  a  real  form,  not  in  a  complex  form.  A  real  and  closed  form 
solution  for  the  Green’s  function  is  very  useful  in  applications. 

A  by-product  of  the  paper  in  [3]  is  the  problem  of  a  composite  wedge 
subjected  to  a  line  force  at  the  apex  of  the  composite  wedge.  The 
composite  wedge  now  consists  of  any  number  of  homogeneous  piezoelec¬ 
tric  wedges  of  different  wedge  angles.  Again,  the  solution  can  be 
expressed  in  a  real  and  closed  form. 

When  the  number  of  the  wedges  in  the  composite  space  or  composite 
wedge  becomes  infinite,  the  material  is  an  angularly  inhomogeneous 
material.  This  means  that  the  material  property  depends  on  the  polar  angle 
0  in  a  cylindrical  coordinate  system  (r,  0,  z ).  The  Stroh  formalism  does  not 
apply  to  an  inhomogeneous  materials,  let  alone  the  special  case  of  angularly 
inhomogeneous  materials.  A  different  approach  must  be  employed.  This  is 
presented  in  [4].  The  solution  is  again  in  an  explicit  real  form.  Real 
materials  may  not  always  be  homogeneous.  Glass  fibers,  for  instance,  can 


6 


•  * 


be  cylindrical  anisotropic.  This  is  a  special  case  of  angularly  inhomoge¬ 
neous  anisotropic  materials. 

The  papers  [3,4]  inspire  the  Russian  scientist.  Professor  V.  I.  Alshits 
and  the  French  scientist.  Professor  H.  O.  K.  Kirchner.  They  pointed  out 
that,  from  a  different  angle,  the  results  in  [3,4]  can  be  re-derived  by  a 
simpler  approach.  In  fact  the  alternate  derivation  also  allows  one  to  include 
piezomagnetic  and  magnetoelectric  properties.  This  results  in  a  joint  paper 
[5],  The  contribution  of  the  principal  investigator  on  this  paper  is  the 
derivation  of  governing  differential  equations  employing  a  dual  coordinate 
system.  The  dependent  variables  are  referred  to  a  rectangular  coordinate 
system  while  the  independent  variables  are  refereed  to  a  cylindrical  coordi¬ 
nate  system.  The  employment  of  a  dual  coordinate  system  follows  from  a 
general  dual  coordinate  system  presented  in  a  book  by  the  principal  investi¬ 
gator  on  "Anisotropic  Elasticity:  Theory  and  Applications"  currently  under 
printing  by  Oxford  University  Press. 

The  last  area  of  applications  under  this  project  is  the  problem  of  a  defect 
in  a  material.  A  defect  may  have  the  form  of  a  void  or  an  inclusion  of 
different  material.  We  consider  the  void  to  have  the  shape  of  an  ellipse, 
and  the  inclusion  to  be  an  rigid  elliptic  inclusion  or  an  inclusion  of  dissimi¬ 
lar  material  [6].  A  crack  is  a  special  case  of  an  elliptic  hole  when  the  minor 
axis  of  the  ellipse  becomes  zero.  Various  problems  are  studied.  For  the 
case  of  an  elliptic  hole,  explicit  real  form  solutions  are  obtained  for  the 
stress  along  the  elliptic  hole  boundary  subjected  to  an  arbitrary  prescribed 
boundary  condition  at  the  elliptic  hole  boundary.  For  an  elliptic  rigid 
inclusion  subjected  to  a  line  force,  a  torque,  and  a  line  charge,  a  real  form 
solution  at  the  interface  between  the  piezoelectric  material  and  the  rigid 
inclusion  is  obtained.  Also  obtained  is  the  general  solution  for  an  piezo- 
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electric  inclusion  of  dissimilar  materials  under  a  uniform  loading  at 
infinity.  The  stresses  along  the  interface  between  the  two  materials  are 
obtained  explicitly.  In  particular,  the  stress  concentration  along  the 
interface  can  be  deduced  directly  from  the  solutions.  These  are  useful 
information  in  applications  in  which  one  is  interested  in  the  stress 
localization  in  the  material  due  to  a  defect,  and  in  the  possible  debonding  of 
a  fiber  in  a  composite. 

In  summary,  we  have  generalized  the  Stroh  formalism  for  anisotropic 
elastic  materials  to  include  the  effects  of  heat  flow  and  piezoelectric 
property.  The  generalized  formalism  allows  us  to  extend  many  solutions 
for  purely  anisotropic  materials  to  materials  subjected  to  heat  flow  or  has 
piezoelectric  effects.  Important  results  are  obtained  for  Green's  functions 
for  the  infinite  space,  composite  space,  and  the  infinite  space  with  an 
elliptic  hole  or  a  rigid  inclusion.  Equally  important  are  the  stress  concen¬ 
tration  around  an  elliptic  hole,  elliptic  rigid  inclusion,  or  an  elliptic  inter¬ 
face  between  two  dissimilar  piezoelectric  materials  subjected  to  various 
loading  on  the  materials.  These  are  very  useful  in  applications,  especially 
in  composite  materials. 
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Abstract.  The  two-dimensional  problem  of  an  elliptic  hole  in  a  solid  of  anisotropic  piezoelectric 
material  is  studied.  The  Stroh  formalism  is  adopted  here.  Real  form  solutions  are  obtained  along 
the  hole  boundary  in  the  case  of  an  arbitrarily  prescribed  vector  field  on  the  hole  surface.  For 
an  elliptic  rigid  inclusion  of  electric  conductor  subjected  to  a  line  force,  a  torque,  and  a  line 
charge,  a  real  form  solution  at  the  interface  is  obtained.  Finally,  general  solutions  for  an  elliptic 
piezoelectric  inclusion  with  uniform  loading  at  infinity  are  investigated. 

1.  Introduction 

In  1958  and  1962,  Stroh  elaborated  the  work  of  Eshelby  et  al.  (1953)  on  two-dimensional 
problems  of  general  anisotropic  elasticity  involving  dislocations,  line  forces,  and  steady  waves. 
This  powerful  and  elegant  approach  was  named  the  Stroh  formalism. 

In  1975,  Barnett  and  Lothe  extended  Stroh’s  1962  paper  to  include  the  piezoelectric  effect 
in  which  an  eight-dimensional  framework  had  been  developed.  Here,  we  consider  the  two- 
dimensional  problem  of  an  elliptic  hole  in  a  solid  of  anisotropic  piezoelectric  material.  Similar 
problems  had  been  studied  by  Pak  (1992)  and  Sosa  (1991).  Although  some  useful  solutions  had 
been  derived  in  these  two  papers,  they  were  both  restricted  to  the  transversely  isotropic  situation. 
In  Pak  s  1992  paper,  special  remote  loading  conditions  were  employed  and  the  concentration 
effect  was  studied.  Likewise,  only  remote  loadings  were  considered  in  Sosa’s  1991  paper. 

Here,  solutions  of  an  arbitrarily  prescribed  loading  on  the  hole  surface  are  derived. 
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Furthermore,  in  the  case  of  an  elliptic  rigid  inclusion  of  electric  conductor  subjected  to  a  line 
force,  a  torque,  and  a  free  line  charge,  real  form  solutions  along  the  elliptic  interface  are  obtained 
which  could  be  used  to  examine  the  concentration  effect.  Finally,  we  investigate  the  situation 
of  an  elliptic  piezoelectric  inclusion  with  uniform  loading  at  infinity. 

In  the  following  basic  solutions  of  the  Stroh  formalism  with  the  piezoelectric  effect  are 
given.  Some  boundary  conditions  are  shown  in  section  2.  In  sections  3  and  4,  a  few  useful 
relations  are  derived.  General  field  solutions  to  the  elliptic  problem  are  obtained  in  section  5 
with  emphasize  on  solutions  along  the  elliptic  boundary.  Such  boundary  solutions  could  be 
employed  to  investigate  the  concentration  effect.  However,  arbitrary  constant  vectors  are 
involved  and  remain  unknown.  They  will  be  determined  in  sections  6,  7,  and  8  in  which 
different  boundary  conditions  are  applied. 

In  a  Cartesian  coordinate  system  (x1 ,  x2,  xs)  the  constitutive  equations  for  piezoelectric 
materials  are  given  by  (Tiersten,  1969) 

°ij  =  CijhnUk,m  +  Di  =  eikmUk,m  ~  H  J  ,k,m  =  1 , 2 ,  3)  G-D 

in  which  repeated  indices  mean  summation  and  a  comma  stands  for  partial  differentiation.  crtJ  is 
the  elastic  stress  and  Di  is  the  electric  displacement.  Coefficients  Cijkm ,  emi- ,  coim  are, 
respectively,  the  elastic  stiffnesses,  piezoelectric  constants,  and  permittivities  with  the  following 
symmetries: 

C ...  =  C...  =  C,  e  =  e  ...  o).  =  o>  ..  (1-2) 

ijkm  jikm  kmij 9  mij  mji  ’  im  mi 

uk  is  the  elastic  displacement  and  <p  is  the  electrostatic  potential.  Cijkm  and  coim  are  positive 
definite  in  the  sense  that 

CijkmUiJUk,n,  >  °»  aimEiEm  >0  (13) 

for  arbitrary  real  nonzero  and  Ei  with 

Et  =  -<P(I.  a-4) 


In  the  absence  of  body  forces  and  free  charges,  the  balance  laws  require 
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%  =  0.  Dut  =  0.  (1.5) 

For  two-dimensional  deformations  in  which  uk  and  <p  depend  on  xI  and  ;c2  only,  a  general 
solution  to  (1.5)  is  given  by 

Uj  =  aj{z)  (J  =  1 ,2,3,4)  (1.6) 

in  which 

z=Xj+px2,  u4  =  <p,  (1.7) 

and  p,  Uj  are  constants  to  be  determined.  In  matrix  notation, 

u  =  a f(z).  (1-8) 


Thus  u,  a  are  four- vectors  and  u  is  called  the  generalized  displacement.  By  defining 


QE 


en  ~Cin 


RE  e21 

>  R  =  T  ,  T  = 

ei2  ~  ^12 


e22  0)22 


where 


(QE)i»  *  cmi’  (*\-CiIk2,  (T %-Cat2,  (eij)w  =  eQm>  d-10) 

we  combine  (1.1),  (1.5),  and  (1.6)  into  one  equation  as 

Q  +p(R  +Rt)  +p2T]a  =  0.  d-11) 

The  4x4  matrices  Q  and  T  are  symmetric  but  not  positive  definite.  However,  they  can  be 
shown  to  be  nonsingular. 

Let  the  generalized  stress  function  vector  $  be  defined  as 

4»  =  bXz),  b  =  (RT+pT)a  =  —  (Q+pR)a,  (1.12) 

P 
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The  second  equality  in  (1.12),  follows  from  (1.11).  Equation  (1.13)  provide  all  components 
of  oi}  and  Dt  except  o33  and  D3  ;  they  can  be  determined  from  (1.1). 

With  the  positive  definiteness  of  Cijkm  and  coim  shown  in  (1.3),  the  eigenvalues  p  of 
(1.11)  are  all  complex  and  consist  of  four  pairs  of  complex  conjugates.  Let 

Pa+4  ~  Pa’  Im{Pa)  >  0  (*  =  1,2,  3, 4),  (!-14) 


a 


d+  4 


cc±  4 


(1.15) 


where  the  overbar  denote  the  complex  conjugates.  The  general  solution  obtained  by  superposing 
eight  solutions  of  (1.8)  and  (1.12),  associated  with  the  eight  eigenvalues  pa  are 

♦  -  (U6) 

in  which  Re  stands  for  the  real  part  and  /ff+4  =Ja  (a  =  1,2, 3, 4)  is  chosen. 


In  most  applications 

fa(za)  =  (a  not  summed)  (1-17) 

is  assumed.  Hence,  equation  (1.16)  reduces  to,  in  matrix  notation, 

u  =  2Re{A<Azm)>q),  <|>  =  2Re{B  </(zJ>  q}  (1.18) 

where  A  and  B  are  4x4  matrices  given  by 

A  =  [a;,  a2,  a2,  a4j,  B  =  [b;,  b2,  b2,  b4],  (1.19) 

and  >  is  the  4x4  diagonal  matrix 

<AO>  =  diag<  f{Zl),  Az2),  Az3),  AzJ  >.  d-20) 


u  =  2Re\YJ  a Ja(za) 


<r=l 


The  elements  of  the  four-vector  q  are  qa  (a  -  1,2, 3, 4).  Notice  that  the  solutions  given  in 
(1.18)  are  in  terms  of  the  arbitrary  function  Aza )  and  the  arbitrary  complex  constant  vector  q. 
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2.  Boundary  Conditions 

Consider  an  arc  or  a  contour  C  described  by 


c(5) :  u = i*/(s)]2+ r  = 


(2.1) 


where  s  is  the  arc-length.  The  unit  tangential  vector  n  and  the  unit  normal  vector  m  are 
given  by 


(2.2) 


respectively.  By  taking  derivative  of  <{>  in  the  direction  of  increasing  s  (with  material  on  the 
RIGHT-hand  side)  and  using  (1.13),  we  obtain 


nT  = 

dxj 

„  j 

dx. 

— ,  0 

,  mT  = 

dx7 

dx. 

—L ,  0 

ds 

ds 

ds  ’ 

ds 

dA  =  , 


d<f>4 


t  {j -1,2,3),  =  D  m  =  D  , 

ds  1  ds  m 


(2.3) 


in  which  tj  is  the  component  of  surface  traction  vector.  Similarly,  one  obtains 

du . 

-T '  =  E  n  =  -E  . 
ds  n 


(2.4) 


If  we  consider  a  dielectric  interface  with  materials  indicated  by  "1"  and  "2",  the  electrical 
conditions  at  the  interface  are 

E,  -n  =  E2  -n,  Dj  m,  +D2  -m2  =  as  (2.5) 

where  n  is  a  unit  vector  tangential  to  the  dielectric  interface,  ntj  is  an  inward  normal  unit 
vector,  and  os  is  the  free  surface  charge  density  along  the  interface.  Without  loss  in  generality, 
we  can  rewrite  (2.5)  as 

<Pi  =  <p2,  D,-m1+D2-m2  =  as.  (2.6) 

If  we  have  an  interface  between  electric  conductor  "1"  and  dielectric  "2",  then  inside  the 


electric  conductor, 
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Dj  =0,  El  =  0. 

In  the  dielectric,  at  the  interface 

E2  n  =  £2n  =  °»  D2  m2  =  <V 


3.  Eight-Dimensional  Formalism 

The  two  equations  in  (1.12)2  can  be  rewritten  as 


f - 

P4 

i 

a 

=  p 

T 

0 

a 

-Q  0! 

b 

R 

I 

b 

Since  T'1  exists,  we  can  reduce  (3.1)  to 

N5  =pl. 


(2.7) 


(2.8) 


(3.1) 


(3.2) 


where 


N,  =  -T->RT,  N2  =  T'1,  N3  =  RT  iRT  -Q  (3.4) 

The  real  4x4  matrices  N2  and  N3  are  symmetric.  Equation  (3.2)  is  a  standard  eigenrelation 
in  the  eight-dimensional  space.  There  are  eight  eigenvalues  pa  (a  =  1,2,...,  8)  and  eight 
associated  eigenvectors  la.  The  eigenvalues  are  the  roots  of  the  determinant 

||N-pI||  =  0.  (3-5) 


N  = 


N,  N2 


n3  N[ 


The  vector  £  in  (3.2)  is  a  right  eigenvector.  The  left  eigenvector  ■q  is  defined  by 

tjtN  =  pr\T,  NTq  =pq,  (3-6) 


and  can  be  shown  to  be 
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T1  = 


(3.7) 


Normalization  of  £a  and  T|^  (which  are  orthogonal  to  each  other)  gives 

nJX  =  dPa  (3'8) 

where  6pa  is  the  Kronecker  delta.  Making  use  of  (1.15),  (1.19),  (3.3)2,  and  (3.7),  equation  (3.8) 
is  written  as 


(3.9) 


This  is  the  orthogonality  relation.  The  two  8x8  matrices  on  the  left  hand  side  of  (3.9)  are  the 
inverses  of  each  other.  Their  product  commutes  so  that 


(3.10) 


BT  AT 

A  A 

1 

I  0 

Bt  At 

B  b 

0  1 

A  A 

BT  AT 

I  0 

B  b 

Bt  At 

0  I 

This  is  the  closure  relation  and  is  equivalent  to 

ABt  +  ABT  =  I  =  BAt  +  BAt,  AAt  +  AAt  =  0  =  BBt  +  BB1 

Hence,  the  three  matrices  S,  H,  L  defined  by 

S  =  i(2ABT  - 1),  H  =  i2AAT,  L  =  -i2BBT 


(3.11) 


(3.12) 


are  real.  The  matrices  H  and  L  are  symmetric  and  nonsingular  (Lothe  and  Barnett,  1976). 
Since  S,  H,  L  are  real,  the  following  relation  exists  (Chung,  1995;  Ting  and  Yan,  1991) 

(3.13) 

which  indicates  that  SH  and  LS  are  anti-symmetric.  It  can  be  shown  that  H^S  and  SL'1 


S  H 

S  H 

I  0 

i - 

H 

H 

xn 

i 

_ i 

0  I 

are  also  anti-symmetric. 

Finally,  we  rewrite  (3.2)  as 


8 


N 

a f'(z) 

a  pf'(z) 

1  i 

b ?(z). 

I 

bpf'iz). 

(3.14) 


Employing  (1.7),,  (1.8),  and  (1.12),  leads  to  a  matrix  differential  equation  for  u  and 


u,i|  _ 

U,2 

4>,2. 

(3.15) 


4.  The  Integral  Formalism 

Let  the  tensor  EUKm  be  defined  by  (Barnett  and  Lothe,  1975;  Kuo  and  Barnett,  1991) 


EuKm  "  Cijkm  UK*  1,2,3), 


^ mij 

eikm 
=  -C0: 


(J  =  1,2,3;  K  =  4), 
(J-4;  K  =  1,2,3), 
(. J  =  K  =  4 ). 


With  n(ca)  and  m(<y)  given  by 

nT(ca)  =  [cos<y,  sinca,  0],  mT(ca)  =  [-sinca,  cosca,  0] 

in  which  a>  is  a  real  parameter  range  from  0  to  2n ,  we  let 

QraC")  =  EUKm  nm^)  ’  RJ^)  =  »|(")  EUKm  ’ 

Tj^o)  =  m£o)Euatmm{6>), 


(4.1) 


(4-2) 


(4.3) 


and 


N \(a>)  =  -T_1(ca)RT(ca),  N2(g>)  =  T'^ca), 

N3(ca)  =  R(o))T\(o)Rt(o))  -Q(o>). 


Lothe  and  Barnett  (1976)  have  shown  that 


7t  TT  71 

S  =  -  JN  x(d)da,  H  =  —J  N 2(a>)da>,  -L  =  —  fN3((o)dco. 

^  n  ^  n  TC 


(4.4) 


(4.5) 


Equations  (4.5)  provide  an  alternate  to  (3.12)  for  the  Barnett-Lothe  tensors  S,  H,  and  L. 
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5.  The  General  Solutions 

In  this  section  the  general  solutions  for  two-dimensional  deformations  along  an  elliptic 
boundary  will  be  derived  (Ting  and  Yan,  1991).  An  ellipse  r  given  by 


j  x[  =  acosifr 
1  x2  =  bsinifr 


(5.1) 


is  shown  in  Fig.l.  Let  n  and  m  be  the  unit  vectors  tangential  and  normal  to  the  elliptic 
boundary,  respectively  and  be  the  angle  between  vector  n  and  the  positive  x1  axis.  Hence, 

nT  =  [cosca,  sinca,  0],  mT  =  [-sinca,  cosca,  0]  (^.2) 


which  is  (4.2).  The  infinitesimal  arc-length  ds  of  the  ellipse  is  given  by 

ds  =  p(i(r )  =  \/tf2sin2  ip  +  b1  cos2  iff . 

From  (2.2),  and  (5.1)  we  see  that 

-a  .  b 

cosca  =  - sm  ilr,  smca  =  - cosi lr 

P(f)  Pit) 


(5.3) 


(5.4) 


when  comparison  is  made  with  (5.2),. 

If  there  is  a  line  force  f  and  a  free  line-charge  density  A.  applied  at  the  origin  (Fig.l), 
by  employing  (2.3),  the  equilibrium  conditions  give 

f  QJjds  =  ^  -  4>/A)  =fj  V- 1.2,3) 

p  B-A 

fD-m ds  =  lim  <f>4(B)  -  <p4{A)  =  -X 

p  B~*A 

in  which  tm  and  D  are  the  surface  traction  and  the  electric  displacement  of  the  medium  along 
the  elliptic  boundary  r,  respectively.  Therefore,  we  have  a  jump  in  <f)  across  the  positive  x; 
axis  if  f  and  X  are  not  equal  to  zero.  Points  A  and  B  are  in  fact  the  same  point  on 
positive  Xj  axis  except  that  when  one  moves  from  A  to  B  counter-clockwise,  the  whole 
ellipse  r  is  transversed. 

Consider  the  transformation 
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=  caCa+daCa-'  (a  =1,2,3, 4), 


(5.6) 


where  ca  and  da  are  complex  constants  and  za=x1  +  pax2.  The  constants  ca  and  da  are 
chosen  such  that  when  (xr  x2)  c  r,  (a  (a  =  1 ,2,3,4)  is  on  a  unit  circle.  That  is, 


Ca  |  =  =  cos  i/f  +  isinf  ( a  =  l,2,3,4 ) 


(5.7) 


when  za  =  acosifr  +  pabsmi/r  and  one  obtains 

a-ipb 


Ca  = 


d*  - 


a  +  ipab 


(5.8) 


Since  a,  b,  and  Im{pa\  are  all  positive  and  non-zero,  it  can  be  shown  that  the  branch 
points  Ca  of  the  transformation  (5.6)  are  located  inside  the  unit  circle  in  the  (a~  plane.  Hence, 
the  branch  points  in  the  (x;,  x2)  plane  are  located  inside  the  ellipse.  In  addition,  the 
transformation  is  one-to-one  outside  the  elliptic  hole. 

In  order  to  satisfy  the  jump  conditions  stated  in  (5.5)  along  the  positive  x{  axis,  the 
arbitrary  function  fiza)  given  in  (1.18)  is  chosen  to  be 

Aza)  =  In  Cm,  (5-9) 


with  za  and  Ca  being  related  by  (5.6).  Also,  by  putting  (Ting,  1986;  1988a,  Hwu  and  Ting, 
1989) 

q  =  ATg0+BTh0  (5-10) 

in  (1.18)  where  gg  and  h0  are  real  constants,  we  obtain  the  first  basic  solution 

u1  =  2ifc(A<taf.>AT}g0  +  2J&{A<lnC>BT}h0  1 
4»'  =  2/?e{B<ln^>AT}g0 +2i?e{B<ln^>BT}h0 

in  which  <  ln^,  >  is  the  diagonal  matrix  of  ln^  with  a  =  1,2,3, 4.  Since  ln^  is  a  multi¬ 
valued  function,  a  cut  along  ifr  =  0  is  introduced  which  makes  u1,  •j*1  single-valued  and  allows 
a  discontinuity  along  the  positive  xl  axis.  As  za=co,  the  elastic  stresses  and  the  electric 
displacements  obtained  from  (5.1 1)2  vanish.  This  is  consistent  with  the  boundary  conditions  at 
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infinity. 

In  order  to  provide  analytical  solutions  outside  the  ellipse, 

Az,)  =  Cl,  q  =  ATgt  +  BTht,  (Jfc-1,2,...)  (5.12) 

are  assumed  in  (1.18)  where  gk,  are  real  constants.  Superimposing  the  solutions 
from  k  =  1  to  oo  leads  to  the  second  basic  solution 

oo  CP 

un  =  2^i?e{A<Ct>AT}g,+2^i?e{A<r:i>BT}hit 

*=1  k= 1 

CO  QQ 

*D  -  2£  Re{B<Ck>Ar\gt  +  2£  fle{B<  f;‘>BT}ht 

^  1  *-  1 

in  which  <(?>  is  the  diagonal  matrix  of  Cak  (a  =  1 ,2,3,4).  Notice  that  both  u11 

and  (I)11  approach  zero  as  za^°°  (or  ^ 

With  (5.7)  it  is  easy  to  see  that 

<la£’.|r>  =  iiff I,  <  C*k\r>  ~  cos(kf)l  -  ism.(ki/r)I.  (5.14) 

Substituting  back  in  (5.11)  gives  the  first  basic  solution  along  the  elliptic  boundary  r  as 

u‘lr  *  *«»,  *‘|r  -  (515) 

=  Hgo  +  Sh0,  (5.16) 

when  using  (3.12).  Similarly,  the  second  basic  solution  along  the  elliptic  boundary  r  is  in  the 
form 

OO  ^ 

09 

u“lr  =  E  [cos (kf)hk  -  sm(ki//)hk],  <frn|r  =  £  \cos(ktfr)gk  -  sm(kifr)gk]  (5-17) 

*=i  *=i  L  J 

in  which 

ht  =  Hg,*Sh„  gt  =  STg,-Lh,.  (5.18) 

Some  useful  relations  between  g,,  ht,  gt,  and  ht  are  given  below  (Chung,  1995;  Ting 
and  Yan,  1991). 


(5.13) 
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K  -  L''(STgt-gt),  K  -  L-‘(gt  +  STg,)  (*  =  0,1,2,...), 
gt=  -H-'(h(tSht),  gt  =  -H-’(Sht-64)  (*-0,1,2,...). 


In  fact,  any  two  of  gk,  h^,  gk,  and  h^.  can  be  written  in  terms  of  the  others. 

In  order  to  investigate  the  concentration  effect,  we  will  derive  the  generalized  stress 
vector  t  and  the  generalized  hoop  stress  vector  tn  along  the  elliptic  boundary.  In  Fig. 2, 
if  n  is  the  arc-length  of  r  measured  in  the  direction  of  n,  then  from  (2.3)  the  generalized 
stress  vector  t  is  defined  as 

m 

=  An]  =  ^  ,n’  (5  —  1) 


or,  using  (5.3),, 


3»l  r 


(5.22) 


Substituting  4*I|Jp  an£l  given  in  (5.15)2  and  (5.17)2  leads 


to 


t1  = 

lm 


1  go,  t“  =  k[sin(lciJf)gk  +  cos(kilr)gk]. 


(5.23) 


Pit) 


Note  that  the  arbitrary  constant  vectors  g^  gk,  and  gk  (k  =  1,2,...)  are  involved. 
Similarly,  in  Fig. 3,  if  the  generalized  hoop  stress  vector  tn  is  defined  by 


*1  -  [CVj.  <‘n)2.  (*.),.  -O.] 


(5.24) 


with  -Dn  =  D  •  (-n),  then  by  letting  m  be  the  arc-length  measured  in  the  direction  of  m,  it 
is  clear  that 


tn  =  =  -4>,;sin6>  +4>2costy 


(5.25) 


where  use  has  been  made  of  (2.3)  and  (5.2)2.  Alternatively,  one  can  express  tn  in  terms  of  u  n 
and  tm  as  (Chung,  1995;  Ting  and  Yan,  1991) 


K  =  N3(<y)u  n  +  N|(6))tm 


(5.26) 
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which  is  a  relation  applies  to  a  general  shape  of  boundary. 

For  the  elliptic  boundary  r  shown  in  Fig.3,  we  have 


K  =  n3(w) 


Mi 


Pit)  df 


NlMi, 


21  2U 

t  +  t 

m  n  n 


in  which 


,n  N-(<y)  “  f  r  A  ) 

tn  =  ‘^^flsin^^h*  +  cos^^hJt  [’ 

nT(«)  A  f  1 

"  -^-E  |^[sin(^ t)  g*  +  cos(k t)  gt]  j 


(5.27) 


(5.28) 


with  the  use  of  (5.15),,  (5.17)!,  and  (5.23).  Again  the  arbitrary  constant  vectors  are  involved. 

In  the  case  of  a  hole  with  free  surface  and  electrically  open  (i.e.,  zero  normal  component 
of  electric  displacement),  (5.26)  then  takes  the  form 

K  =  N3(w)u.n  (5,29) 

and  tm  =  0  (Kuo  and  Barnett,  1991). 

If  we  have  a  rigid  inclusion  of  electric  conductor  inclusion  with  boundary  condition 
E  -n  =  0  given  by  (2.8),,  it  follows  from  (2.4)  that  (5.26)  is  reduced  to 

K  =  N3(w) 

in  which  ur  is  the  rigid  body  motion  of  the  boundary.  Hence, 

ur  =  u0  +  Qe5  x  rr  (5.31) 


+  N[(6>)tn 


(5.30) 


where  uQ  is  a  rigid  body  translation,  Q  is  the  rotation  about  x3  axis  and  rr  is  the  position 
vector  of  a  point  on  the  boundary.  For  an  elliptic  boundary  r,  rr  =  acos^e7  +  bsinifr&2. 
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Thus, 


r/>  =  n>  ur,„  =  Qe3xrr,n  =  Q 


m. 


Substituting  (5.32)-.  into  (5.30)  yields  (Chung,  1995) 


K  =  N[(W)tn 


which  is  also  applicable  to  circular  boundary. 

The  hoop  stress  is  given  by 

a  =  t  • (-n) 

nn  a  V  / 


(5.32) 


(5.33) 


(5.34) 


and  the  two  shear  stresses  are 

anm  =  K  ‘  (~m)>  CTnJ=tn-(-e3)  (5J5) 

in  which  ej  =  [0,  0,  1]. 

Notice  that  our  solutions  are  all  in  terms  of  arbitrary  constant  vectors  gk,  ht,  gk, 
with  k  =  0,1,2,....  When  we  determine  these  constants,  we  have  the  solutions.  In  the  next  three 
sections,  the  arbitrary  constant  vectors  will  be  determined  by  applying  appropriate  boundary 
conditions. 


6. 


An  Elliptic  Hole 

We  consider  an  elliptic  hole  shown  in  Fig.2.  Let 

ea 

ul r  =  “I|r  +  “Ulr  =  +  [cosC**)*^  -  sin(£i/r)h*] 


4>|r  =  tf\r  +  4>U|r  =  f  go  +  S  [c°s(^^)g *  -  sm(kijr)gk] 

jt=i 


(6.1) 


The  right  hand  sides  of  (6.1)  are  given  by  (5.15)  and  (5.17).  Since  u|r  must  be  single-valued, 
it  follows  from  (6.1),  that 
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h0 

and  one  obtains 

ho  =  -S-'Hg,, 


(6.2) 


(6.3) 


when  (5.20)2  is  employed.  The  generalized  stress  vector  along  the  elliptic  hole  boundary  p  is, 
with  (5.23), 


I  ;i  ?n 

t  =  t  + 1  = 

m  mm 


l 


l 


go  -  — —  Y,  UJsin(£^)g^  +  cos (ki/r)gk 


(6-4) 


To  find  tn,  we  first  substitute  (5.19)  into  (6. 1),.  With  (6.2)  and 

SL*1  +L-'St  =  0, 

the  generalized  displacement  vector  along  r  becomes 


(6.5) 


ul r  =  "SL_1E  [cos(^)gfc  ~  sin(kij/)gkj  -  L_1E  [cos(fctfr)gfc  +  sw.(kifr) gfc j .  (6-6) 

Equation  (6.5)  comes  from  the  anti-symmetric  property  of  SL*1.  From  (5.27),  the  generalized 
hoop  stress  vector  then  takes  the  form 


t„  = 


N[(o>)  -N3(<y)SL-1]t, 


N3(o>)L 

P(#> 


m 

-1  oo 


N^SL*1 

p(tf0 


go 


(6.7) 


E  j^sin^i/Og*  -  cos(ktfr)  g^]| 


in  which  (6.4)2  and  (6.6)  are  employed.  Alternatively,  one  obtains 

G3(ca) 


t»  "  G^)tm  + 


P(f) 


i  sT8o  “  E  {*[sin(/:#)g*  “  cos(ki/r)gk] | 


(6.8) 


with 


=  Nj(<u) -N3(ca)SL*1,  G3(w)  =  -N3(<y)L_1 . 


(6.9) 


Since  N3(co)  is  symmetric  so  is  G3(a>)L.  It  can  be  shown  that 
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Gj(6>)L  =  N[(w)L-N3(o>)S 


(6.10) 


is  also  symmetric. 

For  an  arbitrarily  prescribed  boundary  condition  along  r,  we  define  a  four- 

vector 


=  [(TmW)’  (Tm)2(^)’  (Tm AnW] 


[(Of),  DJLt) }  ^ 


in  which  ( i  =1,2,3)  are  arbitrarily  prescribed  traction  components  on  r  while  the 

stress  at  infinity  vanishes.  Dm(ip)  (  =  D(i/r)  -m )  is  an  arbitrarily  prescribed  normal  component 
of  electric  displacement  of  the  medium  along  r  with  the  electric  displacement  vanishes  at 
infinity  also.  Notice  that  Dm(ip)  =  0  refers  to  the  so-called  electrically  opened  situation  (Kuo 
and  Barnett,  1991;  Pak,  1992).  With  the  arbitrarily  prescribed  boundary  conditions  it  is  clear  that 


vw  =  K 


(6.12) 


Employing  (6.4)2  and  the  orthogonality  properties  between  sine  and  cosine,  some  of  the  arbitrary 
constant  vectors  in  terms  of  Tm((/f)  are  determined  as 


2* 

g°  =  p(V)xm{ilr)dxlr 


~r~  ( PWtmWsinikMdf 

if  TT 


0 

It: 


kn 


8/t  =  f  P(^)*m(^)cos(^)<^ 

If  TT  d 


0 

2n 


k  71 


(**  1)  ►. 
(*  *  1  ) 


(6.13) 


Equations  (5.19)2,  (6.2),  and  (6.13),  can  then  be  used  to  find  gQ.  After  that,  h^.  and  h4  are 
obtained  by  employing  (5.19),  and  (5.19)2,  respectively.  In  fact  h0  can  simply  be  computed 
from  (6.3). 

Indeed,  ^  is  related  to  the  resultant  line  force  f  applied  on  r  and  the  resultant  free 
line-charge  density  l  enclosed  by  f.  By  considering  (5.5)  and  (6.12),  the  equilibrium  equation 
becomes 
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f  =  <j)im(t)dn  =  2u 


80 


=  f 


(6.14) 


o  r 

where  n  is  the  arc-length  along  the  elliptic  hole  boundary  r. 

In  the  following  some  special  boundary  conditions  of  xm(f)  are  considered.  We  first 

assume 


p(ifr)D(ip)  =  p  =  constant. 


By  (6.14),  one  has 


DM)  = 


-A, 


p(^r)  2 ti  p(i/f) 


(6.15) 


(6.16) 


Suppose  that  a  uniform  pressure  p  is  applied  along  the  elliptic  boundary  r.  i.e., 

Tm(^)  =  -pm(<»).  (6.17) 


With  this  and  (6.16),,  the  prescribed  generalized  stress  vector  takes  the  form 


= 


p-^-cos^r,  p— —-sin^,  0,  -@— 
Pit)  PW  pit) 


(6.18) 


in  which  (5.2)2  and  (5.4)  are  employed.  By  comparing  with  (6.4)2  the  arbitrary  constant  vectors 
gk  and  gk  are  easily  determined  as 

go  =  ^e4,  gj  =  -pa§2,  g,  =  -p bev  gk  =  g^  =  0  (k  z  2),  (6.19) 


where  the  four-vectors  ^  are  defined  as 


(H=lo;  111 


(6.20) 


Thus,  with  (6.6),  the  generalized  displacement  vector  along  the  elliptic  boundary  is  simply 

u|r  =  SL_1p(x7e2  ~x2ex)  +L'1p(-xie1  +  -jc2e2).  (6.21) 

a  b 


Similarly,  with  (6.8),  (6.12),  and  (6.18),  the  generalized  hoop  stress  vector  is  given  by 
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r  -1 

b 

psincj 

0 

—  COS  03 

a 

-pcosa; 

0 

a  . 

—  sm  a) 
b 

0 

+  G3(6J)ST 

0 

-  G3(£U)p 

P 

p 

0 

p(^)  ( 

.  pit). 

0 

(6.22) 


Likewise,  we  can  consider  a  uniform  in-plane  shear  stress  t  instead  of  a  uniform  pressure  p. 

In  addition  to  the  special  boundary  condition  stated  in  (6.15),  the  traction  boundary 
condition  described  by 


p(ifr)xm(ij/)  =  y  =  constant  = 


2  n 


(6.23) 


is  considered.  Following  the  similar  procedure  given  above,  we  have 

u|r  =  0  =  constant  (6.24) 


which  implies  that  the  elliptic  hole  is  not  distorted  and  the  electrostatic  potential  is  constant  on 
the  hole  surface.  Consequently,  if  the  elliptic  hole  is  filled  with  a  rigid  electric  conductor  and 
subjected  to  a  concentrated  line  force  f  and  a  free  line  charge  density  X  at  the  origin,  the 
generalized  stress  vector  t  along  the  interface  is  simply  given  by  (6.15)  and  (6.23)  (Ting  and 
Yan,  1991).  The  generalized  hoop  stress  vector  is 


— 1— [G,(<i>)  +  G,(w)ST] 

p(«1  1 


Y 

P 


(6.25) 


which  can  be  shown  to  be  consistent  with  what  is  stated  in  (5.33)  when  (6.5)  and  (6.9)  are 
employed. 

In  the  following  the  solutions  for  boundary  conditions  prescribed  as 


dlIml  +  d21m2  +  a3Im 

°  12mi  +  d22m2  +  d32m 
°13mi  +  d23m2  +  d33m 
DjiUj  +  +  D3m3 


3 

3 

3 


D1 


m(w)  = 


a 


Dt 


*«(*>  = 


m(<u) 


(6.26) 
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will  be  derived.  Here,  a  and  f)  are  the  prescribed  uniform  stress  field  and  electric 
displacement  field  along  the  elliptic  hole  boundary  r  within  the  medium,  respectively. 
Since  mT(a>)  =  [-sino>,  cos<u,  0],  we  have  from  (5.4) 


—b  ~ d  ~ 

=  -4- rcos^tj  -  — _  sin^t2,  tx 


p(f) 


°11 

6 12 

°21 

®22 

,  t2  = 

°31 

°32 

A  j 

A 

The  arbitrary  constants  are  determined  by  comparing  with  (6.4)2  which  leads  to 
go  =  0,  =  at2,  gx  =  bit,  gk  =  gfc  =  0  (  k  z  2  ). 

With  the  use  of  (5.1)  and  (6.28),  (6.6)  reduces  to 

u|r  =  -SL  [^^2  -  AC2tj  j  -  L 


b  r  a  7 

— x.t.  +  -xA~ 
a  1  1  b  2  2 


The  generalized  hoop  stress  vector  stated  in  (6.8)  takes  the  form 
tn  =  G^eu)  |cos<yt2  -  sino)^ j  +  G3(o>) 


—  coso>t,  +  — sincuL 
a  1  b  2 


(6.27) 


(6.28) 


(6.29) 


(6.30) 


when  (5.4),  (6,12),  (6.27),,  and  (6.28)  are  employed. 

For  the  problem  of  an  elliptic  hole  subject  to  a  uniform  stress  field  a"  and  a  uniform 
electric  displacement  field  D”  at  infinity  while  the  surface  of  the  hole  is  free  of  traction  and 
electrically  open  (Pak,  1992),  the  solution  may  be  separated  into  two  parts.  The  first  part  is  the 
uniform  solution  in  which  the  stress  and  electric  displacement  are  o"  and  D“  everywhere.  The 
second  part  is  the  "disturbed"  state  due  to  the  presence  of  the  hole.  The  solution  of  the  second 
part  must  satisfy  the  boundary  conditions  that  the  stress  and  electric  displacement  vanish  at 
infinity  while  at  the  hole  surface 


= 


-a 

-D“ 


m(<u). 


(6.31) 


This  is  precisely  the  problem  considered  in  this  section. 
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In  general,  for  an  arbitrarily  prescribed  boundary  condition  Tm(#)  the  series  solutions 
u|r  and  tn  given  above  retain  infinite  terms.  However,  by  introducing  the  conjugate  function 
(Bary.  1964;  Ting  and  Yan,  1991),  one  can  rewrite  the  infinite  series  solutions  in  terms  of 
definite  integrals. 


7.  A  Rigid  Inclusion  of  Electric  Conductor 

In  section  6  the  solutions  for  a  rigid  elliptic  inclusion  of  electric  conductor  in  the  absence 
of  torque  are  studied.  Here,  in  addition  to  a  line  force  f  and  a  free  line-charge  density  X ,  a 
counter-clockwise  torque  Tt3  is  applied.  The  generalized  stress  function  vector  and  generalized 
displacement  vector  along  the  elliptic  boundary  r  are,  respectively, 


Ofl  OP 

<l>lr  =  VZo+Y  \cos(kifr)gk  -sin(^)gj,  «lr  =  Y  [cos(*^)ht  -  sin^tfrjhj.  C7-1) 

*=i  “  J  i=i 


The  equilibrium  conditions  stated  in  (5.5)  are. 


-j*mdn+^  =  0.  f  = 


-X 


(7.2) 


and  t  is  defined  in  (5.21),.  Using  (5. 22)2  and  (7.1),,  (7.2)  reduces  to 


2ti 


-/  ^  pWdV  +  f  =  <|>(0)|r-<|*2rc)|r  +  f  =  -iKgo  +  f-O, 


(7.3) 


J0 

which  is  the  same  result  given  in  (6.14)  so  that  ^  can  be  computed.  Since  the  rigid  inclusion 
does  not  deform  and  the  electrostatic  potential  <p  (  =  u4 )  is  constant  along  the  elliptic  boundary 
jn  by  ignoring  the  constant  components  and  noticing  that  rr  =  acos^re;  +  &sin^e2,  we  have, 
by  using  (5.31), 

u|r  =  Q(acos^re2  -  ftsin^ej).  (7-4) 


Some  of  the  arbitrary  constants  are  determined  by  comparing  with  (7.1)2  which  yields 


hj  =  Qae2,  ht  =  Qbel,  hk  =  hk  =  0  (  lc  z  2). 


(7.5) 
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The  constants  gk  and  gk  are  then  obtained  from  (5.20),  and  (5.20)2,  respectively.  Note  that  the 
angular  rotation  Q  is  still  unknown.  To  determine  it,  we  consider  the  equilibrium  of  moment 
which  gives  (Chung,  1995;  Ting  and  Yan,  1991) 


Q  = 


(7.6) 


where 


U  =  b2H„+a2H£  +  2ab{H~l  S)2I, 
=  c^H-1  +  iH_1S)c  =  ^Zc. 


(7.7) 


In  the  above  cT  =  [-ib,  a,  0,  0]  and  Z  is  the  surface  impedance  tensor  (Lothe  and  Barnett, 
1976;  Kuo  and  Barnett,  1991).  Hence,  one  obtains 


1 


L  +  _Lhi 


tm  = 

m  p(i/f)2n  tzU 


^sin <ue,  +  —  cosoie.  -S(cos£ae,  -sinoie,) 

b  a 


(7.8) 


For  a  circular  rigid  inclusion  p(i/r)  =  a  =  b  =  l,  (7.8)  is  simplified  to 

A  |  A  J1  . 

tm  =  +  — [cosct;e1  +  sino>e2  +  SCsinwej  ~  cos6>e2)].  (7.9) 


In  the  case  of  zero  torque,  one  obtains  from  (7.8) 


A 

_L _ f_ 

p(ilf)  2k 


(7.10) 


This  means  that  p(ij/)im  is  a  constant  which  is  consistent  with  our  observation  in  (6.15)  and 
(6.23).  Furthermore,  for  a  circular  rigid  inclusion,  (7.9)  reduces  to 


V 

=  1 

f 

p » 

2k/ 

-X 

(7.11) 


This  means  that  the  traction  tm  along  the  circular  interface  is  a  constant  in  the  direction 
of  f.  This  phenomena  was  also  observed  in  the  purely  elastic  case  (Ting  and  Yan,  1991). 

To  find  the  generalized  hoop  stress  vector  tn,  (5.33)  is  employed  in  which  tm  is  given 
in  (7.8)  and  (7.9)  for  elliptic  and  circular  rigid  inclusion  of  electric  conductor,  respectively. 
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8.  A  Piezoelectric  Inclusion 

In  this  section  we  consider  an  elliptic  piezoelectric  inclusion  within  a  piezoelectric  matrix. 
If  the  matrix  is  subjected  to  uniform  fields  at  infinity,  with  either  £“  or  a”,  £>“ 
prescribed  such  that  E”  =  o  and  £"  =  0,  then  the  uniform  field  solutions  in  the  absence  of 
an  elliptic  inclusion  can  be  written  as  (Hwu  and  Ting,  1989) 

u~  =  x;e"  +  x2z2,  4>"  =  x&-x2t ",  (8,1) 


in  which 


and 


(8.2) 


(8.3) 


Note  that  (8.1),  is  unique  up  to  a  rigid  body  motion  and  a  constant  electrostatic  potential. 
If  e°°,  E°*  are  given,  o“,  D“  can  be  obtained  by  using  (1.1),  (1.4),  and 

V  7  *  V  * 


2  (uU+uJ.t)- 


If  a"  O'  are  known,  the  following  constitutive  equations  are  employed  (Sosa,  1991): 

eij  ~  ^ijks  ° h;  +  8 kij 7  Ej  -  ~8iksaks  +  Pik^k'  ^  ^ 

The  zero  element  in  implies  that  there  is  no  rotation  of  the  xL  axis  with  respect  to 
the  xs  axis.  By  means  of  superposition  the  field  solutions  in  the  piezoelectric  matrix  are  (Hwu 
and  Ting,  1989) 

u  =  u”  +  2/te{A<<r;1>AT}g1  +22?e{A<(T;1>BT}h1 
4>  =  <|>“  +  2i2e{B<^1>AT}g1 +2^{B<r;1>BT}h1 
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The  first  terms  of  (8.6)  arise  from  the  uniform  fields  applied  at  infinity  without  an  inclusion  and 
are  given  in  (8.1).  The  remaining  terms  are  due  to  the  presence  of  elliptic  inclusion  and  are 
equivalent  to  (5.13)  with  k  =  l .  Thus,  all  the  derivations  obtained  in  previous  sections  can  be 
employed.  Along  the  elliptic  inclusion  boundary  f,  (8.6)  is  reduced  to 

ulr  =  u”|r  +  c°s^hi -sini/r^ 

=  acos^e”  +  ^sin^TE”  +  cosi/rhj  -  sin^-hp 

♦  lr  =  +"lr  +  cos^gi  "  sin^gi  (g  g) 

=  acosi/rt^  -  isin^t”  +  cos^rg1  -  sin^gj, 

in  which  (5.1),  (5.17),  and  (8.1)  are  employed. 

The  solutions  inside  the  piezoelectric  inclusion  are  assumed  uniform  and  have  the  form 

U°  =  xi  6°  +*2e2»  4>°  =  X&~X2t\,  (8-9) 

where 


•-* 

CO 

le°2I  -  Q 

Q 

o  o 

O 

S22 

9  ^  2  — 

2eI3 

2  4 

-e°2 

(8.10) 


*?  -  [°°n.  <4.  <4.  D,°]T  -  -*°2.  «5  -  [<4,  <4,  o°„  D°2 f  =  f,.  (8.11) 

The  constant  Q  in  e°  and  e°  represents  the  rotation  of  the  -^“axis  in  the  inclusion.  Along 
the  elliptic  inclusion  boundary  r,  (8.9)  and  (5.1)  give 


u°|r  =  a  cos  ifre°  +  fcsin^e°,  <|>0|r  =  acos^  -Z>sin^t°. 
The  continuity  condition  states  that 

ulr  =  u°lr’  +1  r  =  4>°l  r 


(8.12) 


(8.13) 
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Employing  (8.7),  (8.8),  and  (8.12)  in  (8.13)  leads  to 


&  E  j  <2  £  j ,  Z?  E  o  h  j  Z>  £2 , 

+  =  at^,  K  +  g,  -  K- 


(B.14) 


In  order  to  solve  h,,  g, ,  e°,  Ej,  t°,  and  t£,  two  more  equations  are  required.  Application  of 
(3.15)  to  the  matrix  and  the  inclusion  yields 


(8.15) 


00 

00 

O 

O 

el 

_ 

e2 

,  N° 

ei 

_ 

e2 

_t" 
ll  _ 

1 

-t? 

With  (5.18)  we  rewrite  (8.14)  in  matrix  form  as 


‘1  ci  ci 

-a  -a 

T  4.O 

>1  to  to 


S  Hlfhj' 

-L  STJ  gl 


(8.16) 


Employing  (3.13)  and  (8.15),  e°,  t^  are  obtained  as 


h  [  s  HI  ei 

—  N°  + 

a  -L  ST  t° 


b  [S  H]  e, 

a  [-L  STj]  £ 


(8.17) 


assuming  that  the  matrix  on  the  left  is  non-singular.  The  constants  hp  gl  and  e^,  t°  are  then 
computed  from  (8.16),  and  (8.15)2,  respectively.  The  rigid  body  rotation  Q  can  easily  be 
determined  from  (8.16),  as 

Q  =  (8.18) 
a 

The  generalized  stress  vector  along  the  elliptic  inclusion  boundary  F  is  given  by 


A 

t  =  d>  = 

m  “n 


asint/ft 2  +  bcosifrt* 


(8.19) 


+  ^^[gisin^  +  8icos^ 


when  (8.8)2  is  used.  Note  that  the  second  part  of  (8.19)  has  the  same  form  as  (5.23)2  with 
k  =  1.  With  (8.8),,  the  generalized  hoop  stress  vector  is 


Observing  that  cos^gj  -  sin^rgj  is  equivalent  to  <j>njr  given  in  (5.17)2  with  k  =  1,  its 
derivative  with  respect  to  m  can  be  computed  as  (Chung,  1995;  Hwu  and  Ting,  1989) 


[cosi/rgj  -sintfrgjj^  =  4>°  =  N3(w)u"  -  N^)^, 


8unL 
=  N,(o>)-  lr 


p(if/)difr 


T,jTr  ^  lr 

pit)difr 


(8.21) 


With  un|r  given  by  (5.17),,  when  k  =  1,  we  have 


t  =  4>  _ 

n  » 


N3(ct»)  r  .  £  i  Nj(w) 

-  hjsm^r  +  hT  cos  f  I  + 


[g1sin^  +  g1  cosier 


(8.22) 


Pit)  L  '  1  '  J  p(ifr) 

Again,  the  second  part  of  (8.22)  is  readily  obtained  by  observing  that  the  result  stated  in  (5.28)2 
can  be  employed  with  k  =  1.  From  (5.4),  (5.25),  and  (8.1),,  <f>“m  is  found  to  be 


4>“m  =  <t>”(-sin<y)  +  <J>^(cosa>)  =  -t£— cos  iff  + 1"— sin  i/r. 

pit)  pit) 


Combining  (8.22)  and  (8.23)  yields,  with  (5.18), 

-1 


(8.23) 


tn  =  pit) 


^bcosif/  -  t”asin^ 


1 


•In, 


pit)  { 


(<y) 


pit)  l 


k 


(«> 


hjSin^r  +  IShj  +  Hg,)cos^r 
g,sin^  +(-Lhj  +STg1)cos^ 


(8.24) 
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Fig.2  Generalized  stress  vector  along  the  elliptic  boundary 


